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Abstract: A fast crack profile reconstitution model in nondestructive testing is developed 
using an arrayed eddy current sensor. The inverse problem is based on an iterative solving 
of the direct problem using genetic algorithms. In the direct problem, assuming a current 
excitation, the incident field produced by all the coils of the arrayed sensor is obtained by 
the translation and superposition of the 2D axisymmetric finite element results obtained for 
one coil; the impedance variation of each coil, due to the crack, is obtained by the 
reciprocity principle involving the dyadic Green's function. For the inverse problem, the 
surface of the crack is subdivided into rectangular cells, and the objective function is 
expressed only in terms of the depth of each cell. The evaluation of the dyadic Green's 
function matrix is made independently of the iterative procedure, making the inversion very 
fast. 

Keywords: arrayed eddy current sensor; superposition principle; ideal crack model; 
reciprocity principle; inverse problem; genetic algorithms 
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1. Introduction 

The use of arrayed eddy current (EC) sensors in Non Destructive Testing (NDT) provides high 
speed inspection and better space resolution by miniaturization of their coils. The arrayed sensors can 
make a measurement of large surfaces without a scan, as illustrated in the Figure 1, which results in a 
gain in time and measurement noise reduction; on the other hand, compared to conventional EC 
sensors, arrayed EC sensors provide more information about the defect characteristics. 

Figure 1. An arrayed eddy current sensor above a piece with a crack. 




There are several configurations of arrayed eddy current sensors [1-3]; when their coils are fed 
separately, the effect of the adjacent coils is negligible; the modeling approach is then the same as for a 
single coil sensors. In this work, we consider an arrayed sensor in which the coils are connected in 
series and fed simultaneously by a current source as shown in Figure 2. The advantages of such 
configuration are: 

The synchronization of the supply and the measurement is not required for the electronic 
component. 

The measurement of the coils impedance is carried throw the voltage measurement. 

The incident electric field on the scan surface is uniform because the coils are connected in 

series, and this is independent of the work piece surface state. 

Figure 2. Impedance matrix measurement. 




Coils matrix nxm 

The investigation is done by the measurement of the impedance variation of each coil. The purpose 
is to determine a crack shape and size using the measurements provided by such a sensor in a real 
time investigation. 

The inversion method proposed is based on the iterative solving of the direct problem; it is thus 
important to have a fast tool to solve the latter. The use of the 3D finite element method would be very 
expensive in terms of memory space and CPU time. On the other hand, the analytic models lack the 
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flexibility to handle complex geometries. In this work, we use the ideal crack model [4-6], generalized 
to arrayed eddy current sensors [7]. In the ideal crack model, the effect of the crack is represented by a 
current dipole layer on its surface, evaluated by an integral equation involving the electric dyadic 
Green's functions and the normal incident electric field on the crack surface. The impedance variation 
of each coil is evaluated using the reciprocity principle. The evaluation of the dyadic Green's function 
matrix is made independently of the iterative procedure of inversion; this makes the inversion to be 
very fast. On the other hand, a fast calculation of the incident field, produced by all the coils of the 
arrayed sensor, on the crack surface is achieved by making a translation and a superposition of the 2D 
axisymmetric finite element results obtained for one coil [7] . 



2. The Modeled System 

Figure 3 describes the modeled system. It is constituted of a (3 x 4) matrix of identical coils situated 
above a conductive plate characterized by a conductivity o and the free space permeability /uo. The 
plate contains an ideal crack of a surface S with an arbitrary shape described in Figure 3. The arrayed 
sensor coils are fed in series by a current source with a time harmonic variation i s (t) = file* 

Table 1 gives the numerical values of the fixed parameters of the system [9]. 



Figure 3. The modeled system. 

db 




Table 1. The Fixed parameter of the modeled system. 



Parameters 



Values 



Frequency: 

Coils: Inner radius, r 2 

Outer radius, r 2 

Height, lb 

Lift-off, If 

Number of turns, N 

Distance between the coils, db 
Plate: Thickness, d 

Conductivity, a 
Crack: Length, L 

Thickness 

Depth 



300 kHz 
0.6 mm 
1.6 mm 
0.8 mm 
0.5 mm 
140 
4 mm 
2 mm 
1 MS/m 
12 mm 
0.2 mm 

Arbitrary shape (Figure 3) 
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Figure 4. Crack shape. 




3. Direct Problem Formulation 



The direct problem is based on the generalization of the ideal crack model to an arrayed eddy 
current sensor [7], which we recall briefly in this section. Firstly, the electric field induced by a single 
coil in the unflawed piece is calculated using the 2D axisymmetric finite element method. The 
electromagnetic problem formulation is given by (1), involving the magnetic vector potential A and the 
current source density Js. The total electric field E T induced by all the coils constituting the arrayed 
sensor is then obtained by (2), making a spatial translation and a superposition of the results obtained 
for the single coil [8]. In (2), Acp is the magnetic vector potential solution of (1) for one coil, nc is the 
number of coils constituting the arrayed sensor, "oxt , oyi? are the center coordinates of the coil k, 
"x,y,z" are the Cartesian coordinates of the computing point, r k is the distance between the computing 
point and the axis of the coil k, and sign(I sk ) indicates the direction of the current in the coil k. 



d 2 8' 



dr dz 



2 rdr 



1 . 

— -J(OMcr 



(1) 



E T X (x, y, z) = -jca£ sign (I sk ) \{r k , z), 



k=l 



E T y (x, y, z) = - si 8 n Usk) * ° Xk \(h . z). 



E T (x,y,z) = 0. 



(2) 



Once the total normal incident field E T N on the surface S of the ideal crack is determined, we 
calculate an equivalent current dipole p normal to this surface S by using the following integral 
equation [9]: 



^(r 0 )-lim{ 7 -^ 0 J s G"(r,rV(^j=0; r 0 



where: 



G nn (r, r') = h.G(r,r').h. 



(3) 



(4) 



In (4), ti is the vector normal to the surface S, and G(r, r') is the electric dyadic Green function 
satisfying Equation (5) and subjected to the same continuity conditions as the electric field. In (5), 
I=xx+yy + zz is the unit tensor and k 2 = jcojU Q a: 
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VxVxG(r,r')-k 2 G(r,r') = IS(r-r') 

The integral equation (3) is solved using the moment method. The crack surface is subdivided into 
(N = n L x n d ) rectangular elements of equal surfaces S e ; the dipole density is considered constant in 
each element. We obtain the following matrix equation: 

{E} = [G]{P}. (6) 

The vectors E and P are of dimension (AO; containing respectively the values of and p for the N 
elements of the crack grid. The matrix G is of dimension (N x AO; its elements are calculated 
as follows: 



g(U)=L G" n (r i ,r j )dS , (i = 1,..N, j = 1,..N). 

J if 3 



(7) 



Using the reciprocity principle, the impedance variation of a coil k of the arrayed sensor is given by 
the following equation: 



AZ*>=-^f s Expels 



(8) 



In (6) E ( n k) and p are scalars representing, respectively, the part of the normal electric field induced 
by the coil k on the surface S, and the normal current dipole solution of (3). The discrete form of (8) is 
given by: 

S N 

AZ (t) =--f-£E(i)P(i). (9) 



4. Inverse Problem 

4.1. Reference data 

The reference data for the inversion are obtained by a 3D finite element computation code 
developed in our laboratory. The computation code is based on the AV-A formulation [10] associated 
to the Gmsh meshing software [11]. We obtained the following impedance variation matrix, 
representing the impedances variations of the (3 x 4) matrix of coils constituting the arrayed 
EC sensor: 



\AZ*\ = 



- 0.0008 - 0.0025i 0.0008 + 0.0008i 0.0009 + 0.0040i - 0.0020 - 0.0005i 

- 0.0464 - 0.0125i - 0. 1204 + 0.0755i - 0. 1783 + 0. 1374i - 0.0835 - 0.0002i 

- 0.0008 - 0.0024i 0.0007 + 0.0008i 0.0008 + 0.004H - 0.0019 - 0.0005i 



(Q) 



Figures 5a and 5b represent the 3D finite element modeled geometry and the 3d plot of \AZ*\ 
respectively; the latter gives an overview of the crack profile. 
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Figure 5. (a) The 3D finite element modeled geometry, (b) The obtained impedances variations. 




(a) (b) 
4.2. Inversion procedure 

The detection of the crack is observed through the variation of the impedance matrix. In the initial 
step, we don't know the exact position and orientation of the crack under the arrayed sensors. The 
adjustment of the position of the latter by looking for the maximum variation of the matrix impedance 
is necessary with the aim of getting the crack in the middle and on the main axis of the arrayed sensor. 
This manual operation makes the inverse problem easier and reduces it to the determination of the 
crack profile. It is assumed that the crack is embedded in a known rectangular area of dimensions Lxd. 
This rectangle is subdivided into N=n L xnj rectangular cells. The crack profile is described by a vector 
q containing n L integer numbers varying between 0 and nj. An example of an arbitrary crack shape 
representation using these discrete values is given in Figure 6. The objective function is expressed 
as follows: 

Sl^tte,-)-* 2 *!,- Qi eN > 0< qi <n d ,i = \...n L ( 1Q ) 

k=l,nc 



Figure 6. Example of a crack shape defined by the discrete values qi. 

qi q„L 









i- 



























The norm used here is the absolute value which takes less computation time than the square root 
norm. For a better consideration of the real part in the minimization of the objective function, we 
separate it from the imaginary part in the impedance variation as follows: 



AZ t (*.)-AZ 4 



Re(AZ t (g,.)-AZ;) 

Re K) 



Im(AZ t (g.)-AZ;) 
MO 



(ID 
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We use the genetic algorithm for the minimization of the objective function (11). Genetic 
algorithms have been widely used, associated to the finite element method for the optimization of 
electromagnetic devices [12,13]. It is based on the principle of natural selection. A set of potential 
solutions (a population) is obtained against fitness criteria, and, through iterations, is refined with 
mutation and recombination [14]. 

In the first step, we suppose that the crack occupies all the rectangle surface (L x d); the matrix G, 
as well as the normal incident field E ( n k) , are then calculated once for all the cells of the grid (n L x n d ) 

independently of the inversion procedure which is explained in Figure 7. It is based on an iterative 
solution of only (6) and (9). In each iteration, the elements q t of the vector q are obtained by the genetic 
algorithm; the matrix G, as well as the vector E ( n k) , are then actualized by eliminating the rows and the 

columns corresponding to the cells which do not belong to the new crack surface; once the actualized 
matrix G(qO is inverted, the vector P is then recalculated by (7) and AZ k is reevaluated by (9). Finally, 
an optimal vector q is obtained; the values of its elements determine the discrete shape of the crack. 
This procedure is very fast since the q t are integer values belonging to a short interval [0 n d ] . 



Figure 7. The inversion flow chart. 
[E] 



Initial 

GnxN 



{qj 



Objective function sfqj 



G(qd -0<HX 



P(qi)=[G( qi) ] '[E] 



AZ(qJ= 

<P(ai),E> 



Genetic 
Algorithm 



Stop 



■ Best {<7,} 



For the considered example we have chosen the discretization (n L x rid= 8 x 4); the variable q t can 
take four values (1, 2, 3 and 4) corresponding respectively to 25, 50, 75 and 100% of the plate 
thickness. For a fast computation, the variable g, is coded on two Boolean variables (00, 01, 10 
and 11). The objective function in finality depends on 2 x n L = 16 Boolean variables which correspond 
in the software implementation to the dimension of one short integer variable. We used the genetic 
algorithm toolbox of Matlab (gatools function) with the following parameters: 



Table 2. The Fixed Parameters for the Genetic Algorithm. 



Parameters 


Values 


Population : 


64 


Crossover rates (Uniform) : 


0.8 


Mutation rates (Heuristic) : 


0.02 
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In Figure 8 we present the inversion results q (hatched cells) obtained for the considered example. 
Since the inversion method is heuristic, the execution time varies from 10 to 40 seconds on a 3.4 GHz 
Pentium 4 PC. The number of generations varies between 50 and 100 depending on the initial solution. 



Figure 8. Inversion results q = [1 21 1 2 3 3 2]. 




The part of computation done by the genetic algorithms is not time consuming according to their 
simple operation. The most consuming time in this inversion is the evaluation of the objective function 
and the inversion of the reduced matrix G(qO. 

5. Conclusions 

We have presented a fast crack profile reconstitution procedure using arrayed eddy current sensor 
data. The use of the reciprocity and superposition principles allows a fast resolution of the direct 
problem. In the inverse problem, which is based on an iterative solving of the direct one, we adopted a 
coarse approximation of the crack profile which is represented by only a few discrete values; this 
makes the inversion procedure very fast when using genetic algorithms. This method can give a real 
time inspection when implemented in an embedded NDT hardware. 
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